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Abstract 

Taking into account both random impurity distribution and thermal fluctua- 
tions of localized spins, we have performed a model calculation for the carrier 
(hole) state in Gai_ x Mn z As by using the coherent potential approximation 
(CPA). The result reveals that a p-hole in the band tail of Gai_ x Mn x As is 
not like a free carrier but is rather virtually bounded to impurity sites. The 
carrier spin strongly couples to the localized d spins on Mn ions. The hopping 
of the carrier among Mn sites causes the ferromagnetic ordering of the local- 
ized spins through the double-exchange mechanism. The Curie temperature 
obtained by using conventional parameters agrees well with the experimental 
result. 
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Although a considerable amount of experimental data has already been accumulated, 
the origin of the ferromagnetism in III-V-based diluted magnetic semiconductors (DMS's) 
has still not been clarified theoretically. Since the magnetic interaction between Mn ions 
has been shown to be antiferromagnetic in n-type (In,Mn)As or intentionally compensated 
(Ga,Mn)As, the ferromagnetic interaction in III-V-based DMS's is most likely hole-induced.S 
The ferromagnetism may be caused by the effective exchange interaction between the delo- 
calized carrier (or the hole originating from the shallow acceptor) and the localized magnetic 
moment ((/spin) on the Mn ion. Matsukura et a/.i referred to the Ruderman-Kittel-Kasuya- 
Yoshida (RKKY) interaction as the origin of the above ferromagnetism. However, the ap- 
plication of the RKKY interaction to this system, as pointed out by various authorsj^ri 
is somewhat questionable. The RKKY formula is applicable when the exchange energy 
IS (= / x S) is small compared to Fermi energy, but this is not the case here because of 
the low carrier density! Furthermore, the measurement of the magnetic circular dichroism 
(MCD) spectrum of (Ga,Mn)As suggests that the holes play more active role in mediating 
the ferromagnetic exchange than traditional RKKY systems.i To date, in regard to the ori- 
gin of the above ferromagnetism, the double-exchange mechanism,! the double -resonance 
mechanismjl the Zener-model description^ and the self-consistent theory based on the CPA 
and mean-field theory! have been proposed. 

In this letter, we apply the dynamical coherent potential approximation (dynamical 
to the carriers in Gai-^Mn^As described by a model Hamiltonian if 

H ^ ] ^m'nP'rauP"np- ^ ] U m , (1) 
m,n,fi m 

where (a m/ J creates (annihilates) a carrier (or p-hole) with spin // at the m-th site. We 
treat Gai_ x Mn x As as an alloy in which a mole fraction x of Ga ions (denoted by symbol A) 
is replaced at random with Mn ions (denoted by symbol M) . Hence U m is Ea <4nu, a mn or 
Em d mtl a mfJi — IJ2^ a l nf i (T ^ ' Sm«mi/ depending on whether the site m is occupied by a Ga 
or a Mn ion. Here, Ea and Em denote the spin-independent potentials for the Ga and the 
Mn ion, respectively; —la ■ S m represents the p-d exchange interaction between the carrier 
and the localized (d-) spin S m of Mn at the site m. Hereafter, we set Ea = as the origin 
of the energy. A carrier moving in a DMS is subject to disordered potentials which arise 
not only from substitutional disorder but also from thermal fluctuations of the localized 
(d-) spins. In the dynamical CPA,liil the disordered potentials are taken into account in 
terms of the spin-dependent coherent potentials, E-f and Sj, which are determined such that 
the effective scattering of a carrier at an arbitrarily chosen site embedded in the effective 
medium vanishes on average. The condition (the dynamical CPA condition) is given by 

(l-x)t^ + x< iff > av = , (2a) 

hj fl + x < *U 

The multiple scattering of an j (j)-spin carrier by an A and an M ion embedded in the 

Note that in the above expression for t*i (iff), the spin flip processes are properly taken 
into account.0 The thermal average over the fluctuations of the localized spins is taken 
by assuming that they are subject to an effective field h, and denoted by < > av . In this 
work we assume the localized spins to be classical. Assuming a semicircular density of 



1 _ x)t A +x< t M >av = . (2b) 

m | (|)-spin carrier by an A and an M ion embedded in the 
effective medium is represented by the t matrix element t& (tf\) and tH (tff), respectively. 
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states(DOS) with the half-bandwidth A for the host valence band, we solved Eqs (|2a]) and 
(|bD simultaneously, and obtained the DOS per site D^(u) (D±(u)) of the carrier with j (|) 
-spin. 

We assume A = 2eV, IS/ A = —0.4 and Em /A = —0.3 as parameters for Gai-^Mn^As; 
-Em (< 0) is the attractive potential acting on a hole due to a Mn 2+ ion (see later discussion). 
In the diluted impurity limit (x — > 0), the present model gives an impurity (acceptor) level 
at the energy of E a = — 1.057A. With increase in Mn concentration an impurity band forms, 
and for x > 0.02 it merges into the host band which originates from the valence band of the 
host semiconductor. This result is consistent with the observed insulator- metal transition in 
Gai-^Mn^As at x ~ 0.03.0 The result of the carrier DOS for x =0.05 is shown in Fig. 1(a) for 
various values of magnetization < S z > av / S. Note that the system with x = 0.05 is metallic. 
The band tail arising from the magnetic impurity band is strongly affected by magnetization. 
With increase of the magnetization the lower edge of the DOS is slightly shifted to the lower 
energy and at the same time D^(u) is strongly suppressed in the band tail. We show the 
optical carrier spin polarization defined by P(u) = (D^(u) — D^(oj))/(D^(uj) + D^(uj)) in 
Fig. 1(b). In order to examine the exchange coupling of carriers with localized spins, we 
calculated the relative local DOS at Mn sites, R(u) = (Df(uj) + Df(u))/(D^(u) + D^u)), 
and depicted it in Fig. 1(c). The result shows that a carrier at the band tail stays mainly at 
Mn sites (R(oo) ~ 0.6) in spite of small x. It should be noted that R(u) is almost independent 
of < S z > av . The result suggests that the change in P(u) is mainly ascribed to the change 
at the magnetic ion sites. 

In order to examine how strongly the carrier spin couples to the localized spins we 
calculate the carrier spin polarization P(n) = (ni — nj) / (nj, + n^) as a function of the total 
carrier density n. Throughout this study, we assume that the carriers are degenerate. In 
Fig. 2, we show the results of P(n) for x = 0.01 and x = 0.05. When x = 0.01, the impurity 
band is separated from the host band, so that carriers stay mainly at Mn sites as long as 
n ^ x. The carrier spin at a Mn site is expected to follow faithfully the fluctuations of 
the localized spin. This leads to the result that P{n) m< S z > av /S for n < x (= 0.01), 
as shown in Fig.2(a). When x = 0.05, the magnetic impurity band partly merges into 
the host band and the Fermi level (e_p) reaches the bottom of the host band at the carrier 
density n\ ~ 0.035. Accordingly, P(n) ~< S z > av /S for n < n\ as shown in Fig.2(b). 
When n^ni, carriers with energy higher than —A seeps out of Mn site to itinerate among 
nonmagnetic sites. 

We calculate the total energy as a function of < S z > av as 

E{< S z > av ) = r u[D^{u) + D^du • (3) 

J — oo 

Note that E(< S z > av ) is the sum of the kinetic and exchange energy. We show in Fig.3(a) 
the energy difference between the ferromagnetic and the paramagnetic state as a function 
of n. Figure 3(a) reveals that the energy gain is approximately proportional to ( <5, ^ >av ^ at 
a fixed n. 

The free energy of the present model per site is given as 

F(< S z > av ) = E(< S z > av ) - TS, (4) 
where the entropy term due to the localized spins is given by 
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s 

S = x/c^log ^2 exp 

s x =-s 



k B T 



hS z 



) 



— x— < S z > 



av 



(5) 



The effective field h is determined so as to minimize F(< S z > av ) through the condition 
J^F(< S z > av ) = 0. If we can expand F(< S z > av ) in terms of (< S z > av ) 2 , T c is determined 
as the temperature where the coefficient of (< S z > av ) 2 vanishes. While calculations for 
small magnetization accurate enough to estimate the coefficient is difficult, we observed from 
Fig. 3(a) that E(< S z > av ) — E(Q)] is approximately proportional to (< S z > av ) 2 up to full 
polarization. Therefore we fitted the data of [E(< S z > av ) — E(0)]/A with (< S z > av / S) 2 = 
0.4 and f .0 to the expansion E(< S z > av ) - E(0)]/A = -a(< S z > av / S) 2 + 6(< S z > av 
/S) 4 and estimated a and b in actual calculation. Then T c is given by ksT c = 2a/3x. In 
Fig. 3(b), the result of T c is depicted as a function of n for x = 0.05. It has been reported 
that ferromagnetism with Curie temperature T c = ffOK is realized in Gai-^Mn^As with 
x = 0.053 when n is 30% of the nominal concentration of Mn.0 The present result is in 
excellent agreement with the experimental result. 

In order to obtain further insight into the mechanism of the carrier-induced ferromag- 
netism in DMS's, we consider the case that the exchange interaction is so strong that the 
carrier spins keep almost antiparallel to the fluctuating localized spins. We show the result 
for IS/ A = —1.0 and Em/ A = in Fig. 4 for illustration. The magnetic impurity band 
is formed about the impurity level (marked 'A') and is separated from the host band. The 
total number of states in the impurity band per site is equal to x irrespective to the value 
of < S z > av / S. When < S z > av = S, all states in the impurity band are down-spin state, 
while when < S z > av = the impurity band is composed of the same number (x/2) of up- 
and down-spin states. The ferromagnetic DOS extends over a wider energy range about 
the impurity level than does the paramagnetic DOS. Consequently the ferromagnetic state 
has lower energy than the paramagnetic state when n is small. This result implies that 
the double-exchange mechanism0 works in the magnetic impurity band to realize ferro- 
magnetism. The energy gain increases first with the increase of n and reaches a maximum 
value at n ~ x/2. Then it decreases gradually and vanishes finally at n ~ x as shown in 
Fig. 3(b). The same mechanism is expected to be operative in the case of IS = — 0.4A 
and Eyi = — 0.3A although the impurity band is not separated in this case. Therefore we 
conclude that the ferromagnetism in Gai-^Mn^As with x = 0.05 is caused by the double- 
exchange mechanism in the band tail. 

Here, we discuss the parameters and the model used in the present work. Since the 
width of the valence band is estimated to be about 4 eV by band calculations,ll!0 we 
have taken A = 2 eV. There are controversial discussions on the magnitude (even on the 
sign) of the exchange constant in III-V-based DMS's.0 Therefore, we have assumed that 
IS/ A = -0.4, which is a conventional value in II-VI-based DMS's.0 Note that \ Ew £ IS \ > 0.5 
is the necessary condition for the existence of an impurity level. Accordingly, the attractive 
potential Em < is needed to reproduce an impurity level, which is confirmed by core- 



reproduces experimental value of the acceptor energy 0.113eV (= 0.057A)B We emphasize 
that the attractive potential strongly assists the occurrence of the ferromagnetism. If we take 
Em/ A = together with IS/ A = —0.4, which corresponds to Mn doped II- VI DMS's, T c is 
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calculated to be at highest 40 K for x = 0.05. This is because the carrier in the II-VI-based 
DMS's spreads over many nonmagnetic sites, which reduces the effective exchange coupling. 
The result explains well why III-V-based DMS's such as (Ga,Mn)As exhibit appreciable 
ferromagnetism while II-VI-DMS's such as (Cd,Mn)Te and (Zn,Mn)Te exhibit only negligible 
ferromagnetism. To the contrary, if we take Em/A- = —0.6 together with IS/ A = —0.4, 
the magnetic impurity band separates from the host band, and T c increases to 140 K for 
x = 0.05. In the latter case, the carrier is semi-localized to a magnetic impurity site by the 
attractive potentials so that the double-exchange mechanism works efficiently 

In the present model, the carrier transfer energy (e mn ) between the Mn site and Ga site is 
assumed to be the same as that between Ga sites. That is, holes in (Ga,Mn)As are assumed 
to have mainly the As 4p character which is originated from the valence-band maximum 
of the host GaAs, rather than Mn 3d character (although we actually treated only carriers 
on cation sites). Therefore, the present model is in contrast to the double-exchange model 
based on the Mn 3d bandJ-3 Results of recent studies using the angle-resolved photoemis- 
sion spectroscopy^ and the band-structure calculation using the local-density-approximation 
with Coulomb interaction (the LDA+U method)0 strongly support our assumption. 

In summary, we have shown that the carriers in the band tail of Gai-^Mn^As have so large 
local-carrier-density on magnetic ion sites that the carrier spins well couples the localized 
spins fluctuating on the Mn ions. The double-exchange mechanism in the magnetic impurity 
band causes the ferromagnetism in the DMS. The attractive local potentials at impurity sites 
strongly assist the occurrence of the ferromagnetism in III-V-based DMS's. 
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FIGURES 



FIG. 1. (a) The Carrier (hole) density of states for x=0.05 are depicted for various val- 
ues of magnetization. The solid and dotted lines represent DAu) and D^(uj), respectively. 
The arrow indicates the Fermi level for n = x (= 0.05), which is about —0.95 A and 
is almost independent of < S z > av . (b) The optical carrier spin polarization defined by 
P{uj) = [Di(u) - D^(u)]/[Di(u) + £>t(w)]. The impurity level {uj/A = -1.057) is also in- 
dicated by a dot (symbol 'A'), (c) The ratio of the local DOS at Mn site to the total DOS, 
R(u) = [Df (w) + Df '(w)]/[D T (w) + D^lo)}, for < S z > av /S = 0.0, 0.5, and 1.0. 

FIG. 2. The carrier spin polarization defined by P(n) = (n^ — n-f)/(nj +n-f) for various degrees 
of magnetization for x = 0.01 (a) and x = 0.05 (b). The arrows indicate n = x and n = n\\ n\ is 
the carrier density where the Fermi level (ef) reaches the bottom of the host band. 

FIG. 3. (a) [E(< S z > ay ) — E(0)]/A is depicted for x = 0.05 as a function of the carrier density 
n, for various values of (< S z > av / S) 2 . (b) The Curie temperature T c calculated for x = 0.05 as 
a function of n. 

FIG. 4. The DOS of the magnetic impurity band in the case of IS/ A = —1.0 and Em/ A = 0. 
The thick, thin and dotted lines show the cases of < S z > av /S = 1.0, 0.5 and 0.0, respectively. 
The magnetic impurity level is dotted with symbol 'A'. 
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